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ABSTRACT 

As part of an effort to associate the unidentified EGRET 7-ray sources with pulsars, 
error boxes for 19 sources were searched using Arecibo at 327-MHz. The sources were 
chosen to be out-of-plane and possibly associated with the Gould belt, a nearby star- 
burst region with an enhanced production rate of core-collapse supernovae. The search 
revealed one new 597-ms pulsar, J2243+1518, within the error box of the EGRET 
source 3EG J2243+1509. The spin-down energy loss rate of the new pulsar is not 
nearly sufficient to power the 7-ray source and so the pulsar is very unlikely to be 
associated. Simulations we have carried out show that any pulsars at Gould belt dis- 
tances should have been detected by the survey. This suggests that either the EGRET 
sources associated with the Gould belt are not pulsars, or that the minimum of the pul- 
sar luminosity function is lower than the ~ 1.5 mJy kpc 2 inferred from the population 
of normal pulsars. 
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1 INTRODUCTION 

One long-standing mystery in astrophysics is the origin of 
the unidentified 7-ray sources. The Energetic Gamma Ray 
Experiment Telescope (EGRET) produced a catalogue of 
271 7-ray sources at 100 MeV to 10 GeV energies (Hart- 
man et al. 1999). These included 169 sources not identifi- 
able with either pulsars or blazars, the two known classes 
of 7-ray source. While some promising identifications have 
been suggested since that time, the majority of these sources 
remain unidentified. As many as 40 of these unidentified 
sources have been shown to be statistically associated with 
the so-called "Gould belt", a dense r egion of atomic and 
molecular gas containing young stars (iGehrels et alJl2000t 
iGrenierlEoobT: lHarding fc Zhandl200 J) . The Gould belt is 
~ 100 pc from the Sun at its nearest point and covers a 
region of radius ~ 300 pc at an angle of nearly 20° from 
the Galactic plane. Its nature as a starburst region means 
it produces core-collapse supernovae at an enhanced rate of 
75 to 95 Myr -1 kpc ~ 2 , ~ 4 times the Galactic rate. Since 
core-collapse supernovae are the most likely mechanism for 
pulsar production, this suggests that t he belt should contain 
a sizable population of young pulsars <lGrenierll2000)) . 

With the possible exception of one millisecond pul- 
sar (iKufper et alJ l200d) . all of the 6-8 EGRET pulsars 
from which 7-ray pulsations have been detected are rel- 
atively young and energetic, with characteristic ages less 
than 500 kyr and spin-down energy loss rates greater than 
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3 x 10 34 ergs s _1 . Several convincing associations between 
newly discovered pulsars and unidentified EGRET sources 
have also been made (e.g. Kramer et al. 2003 and Roberts 
2005). These pulsars are also young and energetic. There- 
fore, most previous surveys of EGRET error boxes for pul- 
sar conterparts have concentrated on those sources close to 
the Galactic plane, where younger pulsars are likely to be 
found. 



The combination of the Gould belt's potential to con- 
tain a population of young pulsars and its statistical asso- 
ciation with the unidentified EGRET sources has prompted 
investigations of EGRET error boxes away from the plane. 
One previous search of mid-latitude error boxes resulted in 
no convincing Go uld belt pulsar dete ctions or EGRET pul- 
sar counterparts jRoberts et alJl2004) . In this paper, we re- 
port on a search of 19 EGRET error boxes using Arecibo at 
a radio frequency of 327 MHz. The sources were chosen to be 
out of the plane and possibly associated with the Gould belt. 
This survey resulted in the discovery of PSR J2243+1518, a 
597-ms pulsar with a spin-down age of 84 Myr. In Section|5| 
we summarize our survey observations and data processing. 
The timing of PSR J2243+1518 is described in Section [3 
We discuss the possible association of the new pulsar with 
3EG 2243+1509 and test the hypothesis that pulsars in the 
Gould belt contribute to the EGRET sources in Section 0] 
Finally, in Section|3we summarize our main conclusions and 
briefly review other possible origins for the EGRET pulsars. 
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2 OBSERVATIONS AND ANALYSIS 

In August 2003 the Arecibo radio telescope was used to 
take the data presented here. The high Galactic latitudes 
ensured that the dispersion measure (DM) of any newly dis- 
covered pulsar would be small (less than ~ 200 pc cm -3 
according to the Cordes & Lazio 2002 model for Galactic 
free electron density). This allowed us to observe at low fre- 
quencies, where pulsar emission is generally strongest. The 
newly commissioned 327-MHz receiver 1 system in the Gre- 
gorian dome was used in combination with the Wideband 
Arecibo Pulsar Processor (WAPP), a fast-dump digital cor- 
relator witii_a_p£ogrammable number of lags and time reso- 
lution dPowd et ail2 000l. In this case 512 lags and a sample 
time of 125 (is were used. The Gregorian dome 327-MHz sys- 
tem has a gain of 11 K Jy , a bandwidth of 25 MHz and a 
system temperature of 113 K. Each pointing was integrated 
for 260 seconds, giving an average limiting detectable flux 
density at 327 MHz of 0.4 mjy to pulsars with a period of 



^30 ms, pulse duty cycle of 10% and a DM of 40 cm 3 



For a typical pulsar spectral index of -1.6 llLorimer et alJ 
I1995T) . the equivalent survey sensitivity at 430 MHz, a fre- 
quency more commonly used for pulsar searches at Arecibo, 
is Sinln,430 = 0.3 mjy. 

We selected a sample of 19 unidentified EGRET sources 
from the third EGRET catalogue with Galactic latitudes 
b > 5° and with declinations in the range —2° < S < 
38° (i.e. visible with Arecibo). These sources are listed in 
Table The 1-a error boxes were gridded in an hexagonal 
pattern with beam centres separated by the beam width. 
The low frequency had the additional advantage of having 
a large beam size (14 ' at 327 MHz) and requiring fewer 
pointings to cover the error boxes. In total, 224 individual 
pointings were observed, with the coverage of the 1-a error 
boxes for all but five sources complete. 

The observations were analysed at the Jodrell Bank Ob- 
servatory using a 182-processor Beowulf cluster (COBRA, 
Joshi et al. 2003). A suite of data analysis tools jLorimerl 
l200lT) was used in combination with scripts to keep each 
node continuously processing from the pool of data. Each 
pointing was dedispersed using 392 trial DMs between 
and 491.2 pc cm -3 . The upper DM limit is over twice as 
high as the expected DMs for any pulsars in this region 
of the Galaxy. The resulting time series were then Fourier 
transformed and searched for periodic signals. Incoherent 
summing of the first 2, 4, 8 and 16 harmonics was used 
to increase sensitivity to narrow pulse profiles. Any result- 
ing candidates with a signal-to-noise ratio (S/N) greater 
than 8 were folded for visual inspection. It took ~ 20,000 
CPU hours to complete the analysis of the 224 pointings. 



3 RESULTS AND FOLLOW-UP 
OBSERVATIONS 

In the ~ 9.6 deg 2 of sky covered one promising candidate was 
found with P = 597 ms, DM ~ 40 pc cm -3 and S/N ~ 17. 
The candidate was re-observed and confirmed as a new pul- 
sar, J2243+1518, in May 2004 using the same observing sys- 




Figure 1. The integrated pulse profile of PSR J2243+1518 at 
327 (left) and 430 (right) MHz. These profiles are for ~ 300 and 
~ 200 minutes of data at 327 and 430 MHz, respectively. Both 
profiles show full pulse phase over 256 bins. These profiles are 
freely available as part of the European Pulsar Network database 
( http:/ /www.jb. man. ac.uk/ ~pulsar/Resources/epn I. 



tem with a 20-min integration. In Fig. we present inte- 
grated pulse profiles at 327 and 430 MHz. 

The pulsar was timed over the next 12 months at 
327 MHz using the WAPPs and at 430 MHz using the Penn 
State Pulsar Machine (PSPM), a 128-channel analogue in- 
terbank spectrometer which samples the incoming voltages 
from the telesc ope every 80 fis over a bandwidth of 7.68 MHz 
JCadwelllll997l) . The timing procedure used was identical to 
that described by Champion et al. (2005). In brief, using a 
preliminary ephemeris, the data were folded modulo the pe- 
riod and cross correlated with a high S/N template profile to 
obtain an accurate time of arrival (TOA) for each observa- 
tion. The TO As were then analysed with the TEMPO 2 soft- 
ware package to produce a phase-connected solution (where 
every rotation of the pulsar is accounted for) over the time 
span of the observations. This resulted in the ephemeris 
given in Table [5] 



4 DISCUSSION 

4.1 PSR J2243+1518 and 3EG 2243+1509 

The spin parameters of PSR J2243+1518, listed in Tabled 
place it in the centre of the normal distribution of pulsars in 
the P-P diagram. Its P and P imply a characteristic age of 
84 Myr, a magnetic field strength of B = 2.6 x 10 11 G and an 
energy loss rate of E — 2.1 x 10 31 erg s _1 . If this pulsar were 
responsible for the 7-ray emission detected by EGRET then 
this emission would have to be powered by its spin-down en- 
ergy loss. The implied efficiency for conversion of spin-down 
energy into 7-rays is r\ = Lj/E = 26 00(47r/)(d/3.5 kpc) 2 , 
where / is the beaming fraction (i.e. fraction of solid an- 
gle swept out by the 7-ray beam), a 7-ray photon index of 
-2 is assumed and the DM inferred distance of 3.5 kpc is 
used (see Table |2J . If the pulsar were converting all of its 
spin-down energy into 7-rays, the distance would have to be 
~ 0.07 kpc for it to be responsible for the 7-ray source. It is 
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Table 1. The selected unidentified EGRET sources accessible to Arecibo with |b| > 5°. 
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In total, 224 independent beam positions were observed, corresponding to a total area of 9.6 deg 2 . 
"The number of pointings taken to cover the error ellipse using the Arecibo telescope at 327 MHz. 
^Coverage incomplete due to time constraints. Number of pointings required appears in parentheses. 
c For a 30-ms pulsar with a pulse width of 3 ms, DM of 40 pc cm -3 and minimum S/N of 8. 



Table 2. The measured and derived parameters for PSR J2243+1518 based on 117 TOAs spanning 494 days. 



Measured Parameters 



Right ascension (J2000) (h:m:s) 


22:43:09.768(4) 


Declination (J2000) (°:':") 


15:18:25.11(11) 


Period (ms) 


596.799464458(4) 


Period derivative (xlO -16 ) 


1.125(7) 


Epoch of period (MJD) 


53258 


Dispersion measure (pc cm -3 ) 


39.828(5) 


327-MHz pulse width (10 per cent) (ms) 


30(2) 


327-MHz pulse width (50 per cent) (ms) 


9(2) 


327-MHz pulse width (eq) (ms) 


15(2) 


327-MHz flux density (mjy) 


0.24(7) 


430-MHz pulse width (10 per cent) (ms) 


29(2) 


430-MHz pulse width (50 per cent) (ms) 


18(2) 


430-MHz pulse width (cq) (ms) 


14(2) 


430-MHz flux density (mjy) 


0.16(5) 


RMS residual to fit (^ts) 


282 


Derived Parameters 


Galactic longitude (J2000) (deg) 


82.8 


Galactic latitude (J2000) (deg) 


-37.4 


Distance (kpc) a 


3.5 


Spin-down age (Myr) 


84 


Spin-down energy loss rate (xlO 31 erg s _1 ) 


2.1 


Surface dipole magnetic field strength (xlO 11 G) 


2.6 



Figures in parentheses are 1-cr uncertainties in the least significant digits as calculated using TEMPO. The arrival times from which 
this ephemeris was derived is freely available on-line as part of the European Pulsar Network (EPN) database 
I http: / / www, jb.man.ac.uk/ ~pu lsar/Resources/epn I. 
"Inferred from the free electron density model of Cordes & Lazio 12002). Distance measurements using this model have a statistical 

error of ~25% although individual errors may be larger. 
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therefore highly unlikely that PSR J2243+1518 is associated 
with 3EG 2243+1509. 

The above arguments suggest that PSR J2243+1518 is 
part of the background population of pulsars at interme- 
diate/high Galactic latitudes. Based on the drift-scan sur- 
veys taken using the Arecibo telescope at 430 MHz, with 
a sensitivity of S m i n , drift ~ 1 mj y we could expect to fin d 
such a pulsar every ~ 40 deg 2 (iMcLaughlin et al]l20o4) . 
Assuming a spectral index of —1.6 to derive Smin, 430 and 
an isotropic population the number of pulsars expected by 
chance is (9.6/40) x (Smm, 43o/S dr ift)~ 3/2 = 1.56. It is, there- 
fore, no surprise that we should find a pulsar which is not 
associated with the EGRET sources in this sample of error 
boxes. 

4.2 Testing the Gould Belt Hypothesis 

To test the hypothesis that the EGRET sources associ- 
ated with the Gould belt are powered by pulsars, a Monte 
Carlo simulation was performed. All of the sources in Ta- 
ble 1 within 15° of the line of sight of the Gould belt 
(modelled as a plane inclined at 20° to the Galactic plane 
with the ascending node at b — 285°) were assumed to 
be associated with the Gould belt and the minimum lu- 
minosity of a detectable pulsar (based on our estimated 
Smin) was calculated. A distance of 0.7 kpc was used, as 
this is the maximum distance of any part of th e Gould 
belt llGuillout et alJll998t lHarding fc Zhandl200ll) . A pul- 
sar was then randomly picked from an assumed luminos- 
ity function with a lower bound of Lmm — 1-5 mjy kpc 2 , 
upper bound of L max = 4.6 x 10 3 mjy kp c 2 and a slope 
of d logiV/d logL = -1 jLvne et al.lll998li . These bounds 
were scaled from the luminosity function derived by Lyne 
et al. to 327 MHz using a spectral index of -1.6. The flux 
density corresponding to this luminosity was then calculated 
and compared to S m i n and the number of detectable model 
pulsars was recorded. In order to minimise statistical fluc- 
tuations in these results, this process was repeated 10,000 
times and the results averaged. 

Of the 19 sources we observed, 12 are within 15° of 
the line of sight of the Gould belt. The simulations al- 
ways detected a pulsar counterpart to each of these sources; 
at a distance of only 0.7 kpc even a pulsar with a lumi- 
nosity of L m i n is det e ctable . Using the beaming model of 
iTauris fc Manchester! il998f) . ~ 2/3 of these pulsars will be 
beamed toward us, and given the coverage of the 1-a error 
boxes, 5.4 pulsars would be expected. Since no pulsars were 
found, it suggests that either the EGRET sources associated 
with the Gould belt are not pulsars or that the L m i n for this 
population of pulsars is lower than L mln = 1.5 mjy kpc 2 
measured in the global population of normal pulsars. For 
example, an L m i n < 0.02 mjy kpc 2 is required before there 
is greater than a 50% chance of not detecting any pulsars. 



5 CONCLUSION 

Our search of 19 out-of-plane EGRET sources has resulted 
in the discovery of one new pulsar, J2243+1518. This new 
pulsar is very unlikely to be the source of the 7-ray emis- 
sion due to the low inferred spin-down energy loss rate. Our 
simulations show that, given the sensitivity of our search 



and the luminosity function of normal pulsars, it is unlikely 
that these error boxes harbour pulsars at distances that an 
association with the Gould belt would suggest. 

Several authors have suggested other origins for the 
mid-latitude unidentified EGRET sources. A possible link 
between Galactic supernova remnants and 19 of the uniden - 
tified EGRET sources is described bv iTorres et"afl J2005I) . 
The 19 positionally coincident sources are at low Galactic 
latitude but it is possible that the high supernovae pro- 
duction rate of the Gould belt could account for some of 
the mid-latitude sources. Low-mass microquasars are poten- 
tial sources for the m id-latitude long-term variable sources 
jGrenier et al J [20051. However these sources are spatially 
distinct from the more stable sources in association with 
the Gould belt. X-ray and radio imaging of unidentified 
EGRET error boxes has resul ted in the discover y of at 
least six pulsar wind nebulae iRoberts et alJ 120051) . How- 
ever, none of these are associated with Gould belt pulsars. 
Conversel y, the sources could be radio-quiet pulsars like 
Geming a jGonthier et alj|2005l) or off-beam 7-ray pulsars 
lHarding fc Zhand feOOliT 

The best opportunity to determine the nature of the 
unidentified EGRET sources wi ll come with next-gene ration 
7-ray telescopes. The AGILE jPellizzoni et aljEobll) tele- 
scope will be able to detect pulsations from the brightes t 
unidentifi ed EGRET sou rces llMcLaughlin fc Cordesll2004h . 
GLAST jThorsettir2003l) . however, will be able to detect 
pulsations from all unidentifie d EGRET sources, unambieu - 
ously determining their origin llMcLaughlin fc Cordesl2 000 ) . 
Furthermore, its excellent angular resolution will allow deep 
radio searches for 7-ray source counterparts to be carried 
out. 
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